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Abstract. The self-consistent quasiparticle RPA (SCQRPA) is constructed to study the effects 
of fluctuations on pairing properties in nuclei at finite temperature and z-projection M of 
angular momentum. Particle-number projection (PNP) is taken into account within the Lipkin- 
Nogami method. Several issues such as the smoothing of superfiuid-normal phase transition, 
thermally assisted pairing in hot rotating nuclei, extraction of the nuclear pairing gap using 
an improved odd-even mass difference are discussed. A novel approach of embedding the PNP 
SCQRPA eigenvalues in the canonical and microcanonical ensembles is proposed and applied 
to describe the recent empirical thermodynamic quantities for iron, molybdenum, dysprosium, 
and ytterbium isotopes. 



1. Introduction 

Sharp phase transitions such as the superfiuid-normal (SN) or shape ones are prominent features 
of infinite systems such as metal superconductors, ultra-cold gases, liquid helium, etc. They are 
well described by many-body theories such as the BCS, RPA or quasiparticle RPA (QRPA). 
The situation changes in finite small systems such as atomic nuclei, where strong quantal and 
thermal fluctuations strongly or completely smooth out these sharp phase transitions. It is 
well known that the conventional BCS, RPA or QRPA theories fail in a number of cases in 
the description of the ground states as well as excited states of these systems. The reason is 
that strong fluctuations invalidate the assumptions, based on which the main equations of these 
theories have been derived. Amongst these assumptions are the Cooper pairs, which violate the 
particle-number conservation, and the closely related quasiboson-approximation (QBA) used in 
the (Q)RPA, which violates the Pauli principle between the fermion pairs. These assumptions 
cause the BCS and QRPA to break down at a certain critical value G c of the pairing interaction 
parameter G, below which the BCS theory only has a trivial solution with zero pairing gap A = 
0. The same is true in the weak coupling region, where the particle-particle RPA is valid but 
its solution also breaks down at G > G c . Meanwhile, the exact solution of the pairing problem 
exposes no singularity at any G [1]. Similarly, at finite temperature T / 0, the omission of 
quasiparticle-number fluctuations (QNF) within the BCS theory leads to the collapse of the 
pairing gap at the critical temperature T c , corresponding to the temperature of the SN phase 
transition in infinite systems. Meanwhile, the exact eigenvalues of the pairing problem embedded 
in the canonical ensemble (CE) shows a smooth decreasing pairing energy with increasing T due 



to thermal fluctuations incorporated in the CE [2]. In rotating nuclei, strong fluctuations also 
smear out the Mottelson-Valatin effect, according to which the pairing gap, existing at zero 
angular momentum M = 0, would collapse at a certain critical angular momentum M c . This 
situation means that, in order to be reliable, the BCS, RPA, and/or QRPA theories need to be 
corrected to include these effects of fluctuations when applied to nuclei, in particular, the light 
ones. This is done within the framework of the self-consistent QRPA (SCQRPA) presented in 
this work. 



2. Formalism 

We consider the pairing Hamiltonian H = J2k>o e kN± k — GJ^kk' -^fc-^V > where N± k = a± k a±k 



is the particle-number operator, and P k 
,t 



(-Pj)t are the pairing operators. The 

operators aj, and a k are respectively the single-particle creation and destruction operators. This 
Hamiltonian has been diagonalized exactly in [1]. The exact partition function is constructed 
by embedding the exact eigenvalues into the CE as ^Exact(/3) = 2~2s d-s ex P(~ /3e| xact ) , with the 
degeneracy d$ = 2 s , inverse temperature j3 = 1/T, and S = 0, 2, ...N being the total seniority 
of the system. Knowing the partition function Z, one calculates the free energy F, entropy 
S, total energy £, heat capacity C, and pairing gap A as F = —T\nZ(T), S = —dF/dT, 
£ = F + TS, C = d£/dT, and A = [-G(£ - 2£ fe e fc / fc + G£ fc /!)] 1/2 , where f k is the 

single-particle occupation number on the kth. level obtained by averaging the state-dependent 

(s) 

occupation numbers within the CE [2]. 

The SCQRPA theory [3, 4] includes a set of BCS-based equations, corrected by the effects of 
QNF, namely 



A fc = A + 5A k , A = Gj2(Vk'}u k >v k , , 5A k 
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where u k and v k are the Bogoliubov's coefficients, 



(2) 



v k 



E k J 2 V E k 

with the renormalized single-particle energies e' k 

G_ 

i. 



'k-Gv 2 k -\\ 

, & k 



e'k ~ Gvl - \y + Al , 



(3) 
(4) 



(V k ) = l — 2n k , the quasiparticle-pair operators A\ = a\a^_ k , A k = and 5M k = n k (l — n k ) 

is the QNF on kth level. To avoid level-dependent gaps A k , the level- weighted gap A^ = 
J2 k A k /Vt [Q, is the number of levels) is considered in the numerical results. Because of coupling 
to collective vibrations beyond the quasiparticle mean field, the quasiparticle occupation number 
n k is not given by a Fermi-Dirac distribution of free fermions, but is found from the integral 
equation [4] 
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where the mass operator M k (oj) and the quasiparticle damping 7fc(w) are functions of n kl 
SCQRPA eigenvalues u^, SCQRPA X£ and y k amplitudes, SCQRPA phonon occupations 
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Figure 1. Energies of the ground state (a) and first excited states for the N = Q = 10 as 
functions of G at T = M = 0. oj pp rpa = £i(N + 2) — £q(N + 2) with the ppRPA eigenvalues £j. 



numbers as well as Uk and Vk- The SCQRPA submatrices A and B contain the screening 
factors and (A\.Ak') so that the set of SCQRPA equations should be solved self- 

consistently with Eqs. (1), (2) and (5) to simultaneously determine A, chemical potential A, 
rik, ujft, X£ and y?. To eliminate particle- number fluctuations inherent in the BCS theory, the 
Lipkin-Nogami (LN) particle- number projection (PNP) [5] is applied on top of Eqs. (1) and (2). 
The ensuing theory, called the LNSCQRPA theory, has also been extended to include the finite 
z-projection M of angular momentum (noncollective rotation) [6]. The set of obtained equations 
is formally the same except that now, depending on the single-particle spin projections ^fjuik 
with 7 being the angular velocity, one has two types of quasiparticle occupation number, n k , 
so that (T>k) = 1 — rik — n_&. At T = and M = the SCQRPA theory reduces to its zero 
temperature and non-rotating limit, where (T>k) = 1/[1 + 2X^(^/fc) 2 ] [3]- 

3. Numerical results and discussions 

Shown in Fig. 1 are the energies of the ground state (a) and first excited state (b) obtained 
at T = M = within several approximations as well as by exactly diagonalizing the pairing 
Hamiltonian for the schematic model, which consists of O doubly-folded equidistant levels with 
the single-particle energies chosen as = k — (Q + l)/2 MeV. The displayed results are for the 
half-filled case with N = £1 = 10, and plotted as functions of the pairing interaction parameter G. 
It is seen that the LNSCQRPA describes rather well the exact energies of both the ground and 
first excited states without any discontinuity in the region around G c , where all other approaches 
such as the RPA, QRPA, and SCQRPA collapse. 

The level- weighted gap, total energy, and heat capacity obtained for the systems with 
N = O, = 10 and 50 are shown as functions of temperature T in Fig. 2. Beside the predictions 
by the SCQRPA, LNSCQRPA, as well as by their corresponding limits, FTBCS1 and FTLN1, 
where coupling to QRPA is omitted (i.e. nt is described by the Fermi-Dirac distribution for 
free fermions), and the finite-temperature (FT) BCS results, the exact CE results are also 
shown. This figure clearly demonstrates how QNF smooth out the sharp SN phase transition 
in finite systems. The pairing gap never collapses, but decreases monotonously with increasing 
T, whereas the spike at T c in the heat capacity, which serves as a signature of sharp SN phase 
transition within the FTBCS, becomes strongly depleted to a broad bump. 

At finite angular momentum M ^ 0, the FTBCS theory predicts the Mottelson-Valatin effect, 
according to which, the zero-temperature pairing gap decreases with increasing M and collapses 
at M = M c because the angular momentum blocks the levels close to the Fermi surface [Fig. 3 
(a) and 3 (c)] . Thermal effects relax the blocking, opening some levels around the Fermi surface 
for pairing. This leads to the thermally assisted pairing gap (or pairing reentrance) , according to 
which at a certain T = T\ the pairing gap becomes finite even at M > M c [7, 8]. With increasing 
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Figure 2. Level-weighted pairing gap A, total energy E, and heat capacity C, as functions 
of T for N = Q = 10 (a - c) and 50 (d - f) obtained within the FTBCS (dotted), FTBCS1 
(thin solid), FTLN1 (thin dashed), SCQRPA (thick solid), LNSCQRPA (thick dashed). The 
dash-dotted lines for N = 10 are the exact CE results. 
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Figure 3. Level-weighted pairing gaps A for N = Q = 10 as a functions of temperature T at 
various values of M/M c and angular momentum M at various values of T/T c within the FTBCS 
(a, c) and FTBCS1 (b, d) theories. 



T thermal effects again break the pairs so that the gap disappears at T = T2 > T\ [See Fig. 3 (a) 
for M/M c > 1]. In finite systems, the QNF smooth out both the Mottelson-Valatin transition 
and thermal assisted pairing, for instance, for N = Q = 10 with G = 0.5 MeV at T/T c > 1, 
the gap only decreases monotonously with increasing M but never vanishes [See Fig. 3 (d) for 
M/M c > 1], whereas at M/M c > 3, the pairing gap reappears at T > T\ but remains finite with 
further increasing T [See Fig. 3 (b) for M/M c > 3]. 

The odd-even mass difference contains the admixture with the contribution from uncorrelated 
single-particle configurations, which increases with T. Therefore, the simple extensions of 
this formula to obtain the three-point and four-point gaps, in principle, do not hold at finite 
temperature. We propose an improved odd-even mass difference formula at T 7^ 0, namely 



S> = ~[(£(N + 1)> Q + (£(N - 1))J - <£(iV)><?> , (6) 



where (£(N)) a is the total energy of the system with N particles within the grand canonical 
ensemble (GCE) (a =GC) or CE (a =C); = 2£ fc [e k - Gf ( k a) /2]f { k a) with -GEkif^? 
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Figure 4. Pairing gaps extracted from the odd-even mass differences as functions of T for N = 
10 (a,c) and N = 9 (b,d) (SI = 10, G = 0.9 MeV). The thin solid and thick solid lines denote 
the gap A® (iV) (A^ A \N) = [A (3 \N) + A&\N - l)]/2), and the improved gap A (i \(3,N) (i = 
3, 4) from Eq. (6), respectively. The dash-dotted lines are the canonical gaps Aq . 




0.5 1.0 1.5 2.0 2.5 0.2 0.4 0.6 0.8 1.0 1.2 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 
T(MeV) T(MeV) T(MeV) T(MeV) 




5 10 15 20 5 10 15 20 5 10 15 20 5 10 15 20 
E*(MeV) E*(MeV) E*(MeV) E*(MeV) 



Figure 5. CE heat capacities C as functions of T and MCE entropies S as functions of excitation 
energy E* for 56 Fe, 94 Mo, 162 Dy, and 172 Yb. Experimental data are taken from [11]. 



coming from uncorrelated single-particle configurations. Shown in Fig. 4 are the pairing gaps 
A^\j3,N) (i = 3 and 4), obtained for iV = 9 and 10 (Q = 10) by using the simple extension 
of the odd-even mass formula to T ^ as well as the modified gaps A^(f3,N) from Eq. (6), 
and the canonical gaps A^. It is seen in Fig. 4 that the naive extension of the odd-even mass 
formula to T / fails to match the temperature-dependence of the canonical gap AqK The 
gap A( 3 )(/3, N = 9) even turns negative at T > 2.4 MeV, suggesting that such simple extension 
of the odd-even mass difference to finite T is invalid. The modified gap AW(/3) is found in 
much better agreement with the canonical one, whereas the modified four-point gaps A^\f3) 
practically coincide with the canonical gaps. The comparison in Fig. 4 suggests that formula (6) 
is a much better candidate for the experimental gap at T ^ 0, rather than the simple odd-even 
mass difference. 

In order to construct a feasible description for pairing within the CE, the eigenvalues of the 
LNBCS and LNSCQRPA, obtained for each total seniority S at T = 0, are embedded into 
the CE by using the partition function Z 7 (£) = T,s d se~ Pe s (7 = LNBCS, LNSCQRPA). The 
resulting approaches are called the CE-LNBCS and CE-LNSCQRPA, respectively [10]. These 
solutions are also embedded into the microcanonical ensemble (MCE) by using the Boltzmann's 
definition for entropy, S(£) = \nW(£), where W{£) is the number of accessible states within 



the energy interval (£,£ + 6£). The corresponding approaches are called the MCE-LNBCS and 
MCE-LNSCQRPA, respectively [10]. 

The CE heat capacities and MCE entropies for several nuclei are shown in Fig. 5 as functions 
of T and excitation energy E* , respectively. The single-particle energies are calculated within 
the deformed Woods-Saxon potentials. In order to have a consistent comparison with the 
recent experimental data in [11], we carried out calculations by using the CE-LNBCS and 
CE-LNSCQRPA for 56 Fe, where pairing is included within the complete pf + gg/ 2 shell above 
the 40 Ca core. For Mo isotopes, pairing is included in the 22 degenerated single-particle levels 
above the 48 Ca core. For Dy and Yb the same is done on top of the 132 Sn core. It is clear from 
this figure that the CE-LNSCQRPA results agree quite well with the experimental data [11], 
which are also deducted from the CE. The MCE entropies, obtained within the MCE-LNBCS 
and MCE-LNSCQRPA using 5£ =1 MeV, are plotted versus the experimental data. It is seen 
that the MCE-LNSCQRPA entropy not only offers the best fit to the experimental data but 
also predicts the results up to higher E* > 10 MeV. 

4. Conclusions 

The proposed LNSCQRPA theory can describe without discontinuity the pairing properties of 
hot noncollectively rotating nuclei at any values of pairing interaction parameter G, temperature 
T, and angular momentum M. In the limit of zero temperature and zero angular momentum, it 
offers the best fits to the exact solutions in the weak coupling region with large particle numbers 
for the energy of the first excited state, whereas the SCQRPA reproduces well the exact one in 
the strong coupling region. In the limit of very large G all the approximations predict nearly 
the same value as that of the exact one. Under the effect of QNF, the paring gaps obtained 
at different values M of angular momentum decrease monotonously as T increases, and do not 
collapse even at hight T. The effect of thermally assisted pairing (pairing reentrance) shows up 
in such a way that the pairing gap reappears at a given T\ > and remains finite at T > T±, 
in qualitative agreement with the results of Ref . [9] . We suggest a novel formula to extract the 
pairing gap at T ^ from the difference of total energies of even and odd systems, where the 
contribution of uncorrelated single-particle motion is subtracted. Its prediction is found in much 
better agreement with the canonical gap than the simple extension of the odd-even mass formula 
to T ^ 0. Finally, we embedded the solutions of the LNBCS and LNSCQRPA into the CE and 
MCE, and found that the CE-LNSCQRPA predictions are in quite good agreements with the 
exact results as well as the recent experimental data. The present approach can also describe 
simultaneously and self-consistently the experimentally extracted total energy, heat capacity, 
and entropy within both CE and MCE treatments. It is simple and feasible for the application 
to larger finite systems, where the exact matrix diagonalization and/or solving the Richardson 
equation are impracticable to find all eigenvalues, whereas other methods, such as the quantum 
Monte-Carlo calculations, are time consuming. 
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